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ABSTRACT 

By obtaining photometric redshift information, tomography allows us to cross-correlate galaxy ellipticities in 
different source redshift bins. The cross-correlation is non-vanishing because the different bins share much of 
the foreground mass distribution from which, over Gpc scales, the lensing signal is built. If the redshift bins 
are thick enough however, the cross-correlations are insensitive to contamination from the intrinsic alignments 
of galaxies since these fall off rapidly on scales larger than a few tens of Mpc. We forecast how lensing to- 
mography using only the cross-power spectra can constrain cosmological parameters compared to tomography 
including the auto-spectra. It is shown that the parameter errors are degraded by only O(10%) for 5 or more 
source redshift bins. Thus, the cross-power spectrum tomography can be a simple, model-independent means 
of reducing the intrinsic alignment contamination while retaining most of the constraints on cosmology. 
Subject headings: cosmology: theory - gravitational lensing - large-scale structure of universe 
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1. INTRODUCTION 

As our cosmological knowledge and techniques mature, 
it becomes increasingly important to strive for observations 
and methods which are, as much as possible, immune to 
known systematic effects as well as being statistically pow- 
erful. Nowhere is this more true than in the field of weak 
gravitational lensing. Since its detection by several groups 
three years ago (Bacon, Refregier & Ellis 2000; Kaiser, Wil- 
son & Luppino 2000; Van Waerbeke et al. 2000; Wittman et 
al. 2000), weak gravitational lensing by large-scale structure 
has become a well established technique, used already to set 
constraints on the mass density (il m o) and the fluctuation am- 
plitude (erg) (e.g., see Van Waerbeke & Mellier 2003 for the 
current status) and touted for its potential to constrain clus- 
ter scaling relations (Huterer & White 2002) and dark energy 
(Benabed & Bernardeau 2001; Huterer 2002; Hu 2002a,b; 
Heavens 2003; Abazajian & Dodelson 2003; Refregier et al. 
2003; Jain & Taylor 2003; Bernstein & Jain 2003; Takada & 
Jain 2003, hereafter TJ03). 

There are many sources of systematic errors which can af- 
fect lensing measurements. One source in particular, intrin- 
sic ellipticity alignments of source galaxies (Heavens, Re- 
fregier & Heymans 2000; Catelan, Kamionkowski & Bland- 
ford 2001; Crittenden et al. 2001; Pen, Lee & Seljak 2001; 
Mackey, White & Kamionkowski 2002; Jing 2002) cannot be 
mitigated by improved instrumentation or reliably predicted 
by theory. While likely not important for the current genera- 
tion of experiments (Heymans et al. 2003), it could prove to 
be a limiting uncertainty in more ambitious ongoing and fu- 
ture surveys such as the CFHT Legacy Survey, the Deep Lens 
Survey, DML/LSST, Pan-STARRS and SNAP. 

In this Letter, we discuss a simple, largely (astrophysical) 
model independent, technique to reduce the susceptibility of 
weak lensing measurements to intrinsic alignments of galax- 
ies. This technique assumes only that the intrinsic alignments 
fall off on scales above a few tens of Mpc while the lens- 
ing signal builds up over Gpc scales. We emphasize the ad- 
vantages of splitting the galaxy distribution into multiple red- 
shift bins and considering only the cross-power spectra. These 
have no shot noise bias, contain most of the cosmological in- 
formation and should be totally insensitive to intrinsic align- 



ments. While not as powerful as model dependent methods 
(Heymans & Heavens 2003; King & Schneider 2002, 2003; 
Heymans et al. 2003) these latter require modeling of intrinsic 
alignments, which involves many uncertain aspects of galaxy 
formation. 

2. TOMOGRAPHY OF LENSING CROSS-POWER SPECTRUM 

2.1. Weak lensing field 

All future surveys will provide photometric redshift infor- 
mation on source galaxies. This additional information allows 
us to subdivide the galaxies into redshift bins which is crucial 
if these surveys are to be used to constrain the evolution of 
cosmological parameters. We shall assume throughout that 
we have the ability to divide our source population into red- 
shift bins, commenting on some of the issues in the conclu- 
sions. 

In the context of cosmological gravitational lensing and 
assuming the Born approximation (Blandford et al. 1991; 
Miralda-Escude 1991; Kaiser 1992, and see Jain, Seljak & 
White 2000 and Vale & White 2003 for tests), the lensing con- 
vergence field can be expressed as a weighted projection of 
the 3D density fluctuation field between observer and source 
(also see Bartelmann & Schneider 2001; Mellier 1999 for re- 
views): 



k(6)= / d X W(x)5[x,X0], 



(1) 



where 9 is the angular position on the sky, x lS the comov- 
ing distance, and xh is the distance to the horizon. We shall 
assume throughout a spatially flat universe. 

The lensing weight function, Wq, for galaxies of subsample 
; sitting in a redshift bin i can be written 



W X 



dXsPs(z) 



dz Xs"X 



X < XM 



(2) 



W( i )(x)=< i nia(x)J Xi ~™" r ° s ~" dXs Xs 

.0, X>Xm- 
where a(x) is the scale factor, p z is the redshift distribution of 
source galaxies and Wq = (3/2)f2 m r>ffg, with Hq the Hubble 
constant and fi m o the present day value of the matter density 
in units of the critical density. Following TJ03, we assume 
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with zo = 0.5 and the average number density of galaxies per 
steradian, n g , which peaks at 2zo = 1 and has median redshift 
Zmed = 1.5. The quantity ;T, is the average number density of 
galaxies in a redshift bin i, defined to lie between the comov- 
ing distances Xi and Xi+i- «< = SxT dXs P^lt > and we 
have assumed sharp subdivisions of the galaxy redshift distri- 
bution for simplicity. 

It is worth noting the dependence on the cosmology in these 
expressions. The overall shear amplitude is sensitive to VL m Q 
and as, which explains the status of cosmological constraints 
derived from current lensing surveys. Tomography allows us 
to recover the redshift evolution of the efficiency, W(j), and 
the mass clustering and the evolution of \ f° r non power- 
law power spectra. Each of these is sensitive to the equa- 
tion of state of dark energy, however the constraints are de- 
termined mainly by the dependence of the lensing efficiency 
(see e.g. Figure 2 in Abazajian & Dodelson 2003). 

2.2. The power spectrum and its covariance 

Using the flat sky approximation and the Limber approxi- 
mation (Limber 1954; Kaiser 1992), the angular power spec- 
trum between the lensing fields of redshift bins i and j is 

Q ; y)(0 = J o XH dxW (i ixW(j)(x)x~ 2 P S (k=^;xj, (4) 

where Pg{k) is the 3D mass power spectrum. The Limber ap- 
proximation holds well over the angular scales we consider: 
50 < / < 3000 (Jain et al. 2000; White & Hu 2000; Vale & 
White 2003). For I > 100 the major contribution to C^(l) 
comes from non-linear clustering (Jain & Seljak 1997; also 
see Figure 2 in TJ03) and we employ the fitting formula of 
Smith et al. (2003) for Ps(k), assuming that it can be applied 
to dark energy cosmologies (e.g. White & Vale 2003). We 
note in passing that the issue of accurate power spectra for 
general dark energy cosmologies still needs to be addressed 
carefully (see Huterer 2002 for related discussion). 

Assuming that the intrinsic ellipticity distribution is uncor- 
rected between different galaxies, the observed power spec- 
trum between redshift bins i and j can be written (Kaiser 
1992, 1998) 

2 

q$(l) = C (iJ) (l) + 6 ij ^, (5) 

where a e denotes the rms of the intrinsic ellipticities and the 
Kronecker delta function accounts for the fact that the cross- 
power measurement is not biased by shot noise. In addition, 
the cross-power is totally insensitive to intrinsic ellipticity 
alignments, if the bins are much larger than 10 Mpc (Heymans 
& Heavens 2003; King & Schneider 2002, 2003; Heymans et 
al. 2003). 

Assuming Gaussian statistics, the covariance between the 
power spectra C(y)(Z) and Q m „)(/') is 

COV = (2/+l)AZ/ sky [ C « C W + C (°») C 0>»)] ' (6) 

where f^ y is the fraction of sky covered, Al is the bin width 
and we have suppressed the argument, I, of C obs . If we re- 
strict our analysis to angular scales I < 3000 the statistical 
properties of the lensing fields are close to Gaussian (White 
& Hu 2000; Cooray & Hu 2001). Further the shot-noise wiU 
provide the dominant contribution to the covariance on the 
small scales where tomography derives most of its cosmo- 
logical constraints, strengthening the case for our Gaussian 



error assumption. It is important to note that even if we con- 
sider only the cross-power spectra, with i ^ j and m^n, the 
shot noise contributes to the covariance when i = m and so 
on. Hence, the shot noise induces the statistical errors in the 
cross-power spectrum measurement, which in turn propagate 
into cosmological parameters. 

Table [2 compares parameter forecasts from tomography 
with and without the auto-spectra. We consider angular scales 
50 < / < 3000, and assume / s k y = 0.1, n g = 100 arcmin" 2 and 
er e = 0.4. Note that all the errors scale with the sky cover- 

-1/2 

age as / sky . Given tomography with n s redshift bins, each 
bin is chosen so as to have an equal number density of galax- 
ies for the redshift distribution (yj; e.g., < zi < 1, 1 < Zi < 
1.7, 1.7 < Z3 for three redshift bins. For all the cases we have 
considered the redshift bins are sufficiently thick, Az ^ 0.2, 
that contamination from intrinsic alignments will be signifi- 
cantly reduced. 

2.3. Fisher matrix analysis and the fiducial model 

Following TJ03, we use the Fisher matrix formalism to ex- 
amine how lensing tomography can constrain cosmological 
parameters. Assuming the likelihood function for the lens- 
ing power spectrum to be Gaussian, the Fisher matrix can be 
expressed as 



dC m (l) 

dp a 



[COV]" 



dC( mn )(l) 
dp/3 



(7) 



where the inverse covariance matrix is denoted by [Cov] -1 
with Cov given by equation (|6) and the p a denote the cos- 
mological parameters. To use only the cross-power spectra 
we simply omit ; = j and m = n 'm the summation of equation 
Q. For n s redshift bins, we can use n s (n s — 1)/2 cross-spectra. 
Hence, to extract redshift evolution of the lensing observables 
requires at least 3 redshift bins. 

The parameter forecasts derived are sensitive to the parame- 
ter space used and to whether constraint on a given parameter 
is obtained by marginalizing over other parameter uncertain- 
ties. We use seven parameters to which the lensing observ- 
ables are sensitive within the cold dark matter (CDM) model: 
^de, wq, w a , C8, «.s Qbh 2 and h, where f^e and are the 
density parameters of dark energy and baryons, n s is the spec- 
tral index of primordial scalar perturbations, h is the Hub- 
ble parameter, and a% is the rms mass fluctuation in a sphere 
of radius 8/! _1 Mpc. We use a simple parameterization of the 
equation of state of dark energy: w(a) = wo+w a (l -a) (Turner 
& White 1997; Linder 2003) and assume a spatially flat uni- 
verse. For the input linear mass power spectrum, we employ 
the BBKS transfer function (Bardeen et al. 1986) with the 
shape parameter given by Sugiyama (1995). 

For the fiducial model we take ACDM model with = 
0.65, w = -1, w a = 0, o- 8 = 0.9, n s = 1, ti b = 0.05 and h = 0.72, 
which is consistent with the recent WMAP results (Spergel et 
al. 2003). We assume the priors o-(lnflbh 2 ) = 0.010, cr{n s ) = 
0.008 and a(h) = 0.13, expected from the Planck mission (see 
Table 2 in Eisenstein et al. 1999). Assuming Gaussian priors, 
we add the diagonal component (fj,rior)a/3 = <5a/3c(Pa)~ 2 to the 
lensing Fisher matrix. 

3. RESULTS 

The table shows that, for three redshift bins, cross-power 
spectrum tomography more than doubles the parameter er- 
rors from tomography including the auto-spectra, because we 
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FIG. 1 . — Projected 68% CL constraints in the parameter space of Q^s, Wfj> and w a for 5 bin tomography with and without the auto-power spectra. 



Parameter Estimation (n g = 100 arcmin 2 , = 0.1, cr e = 0.4) 
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Table 1 

Summary of parameter constraints from tomography of only using the cross-power spectra (Cross-PS) or 
tomography including the auto-spectra (ps). the errors include marginalization over the other parameters. 



can use only 3 of the 6 spectra. However, adding even one 
redshift bin drastically improves the constraints. As a result, 
cross-power spectrum tomography for n s > 5 recovers most 
of cosmological information contained in the auto- and cross- 
power spectra: degradation in the 1-cr errors are as small as 
15%. FigureQshows the constraint ellipses in the param- 
eter space of fide, wo an d w a for five redshift bins. Note that 
the ellipses correspond to 68% confidence level (A% 2 = 2.3), 
including marginalization over the other parameters. In this 
case we can use 10 cross-power spectra for the tomography. 
It is clear that degeneracy directions in the parameter space 
are almost identical between the two tomography methods, 
and the ellipse areas differ by only ~ 30%. 

TJ03 showed that using both power spectrum and bispec- 
trum tomography provides improvements in parameter con- 
straints of a factor of 3 over power spectrum tomography 
alone. This is because significant additional information is 
contained in the non-Gaussian nature of the large scale struc- 
ture inducing the weak lensing. While bispectra can be 
constructed from n s redshift bins for tomography, only the 
n s auto-bispectra are contaminated from intrinsic alignments, 
allowing us to use the other n s (n%— 1) bispectra for cross- 
bispectrum tomography. The column labeled with "+ Bisp" in 
Table ^compares the results when we use all the power spec- 
tra and bispectra or use only the cross-spectra for five redshift 
bins. As expected, cross-bispectrum tomography loses little 
cosmological information and significantly improves parame- 
ter errors derived from the cross-power spectrum tomography. 

4. CONCLUSION AND DISCUSSION 



We have proposed lensing tomography that only uses the 
cross-spectra constructed from source galaxies in different 
redshift bins. Cross-spectrum tomography provides a simple, 
robust and model independent means of reducing systemat- 
ics from intrinsic alignments which retains most of the cos- 
mological information. For more than 4 redshift bins, cross- 
spectrum tomography yields errors on cosmological parame- 
ters only ,$ 15% larger than errors including the auto-spectra 
(see TableQ. The situation is further improved by combining 
with cross-bispectrum tomography. 

As future experiments concentrate on subdividing the 
source galaxies into multiple redshift bins to better reconstruct 
the lensing effect as a function of distance, it becomes more 
important to guard against contamination by intrinsic align- 
ments. Intrinsic alignments are difficult to predict, with cal- 
culations differing by an order of magnitude at present. How- 
ever, simply neglecting intrinsic alignments could bias esti- 
mates of cosmological parameters from future lensing surveys 
(Heymans et al. 2003). Our suggestion is intermediate be- 
tween detailed modeling of the effect and simple neglect. We 
believe it is worth exploring a 3D mass reconstruction tech- 
nique using only the lensing effects on galaxies in different 
redshift bins even if this is only used as a 'sanity check' dur- 
ing the analysis. 

Finally we remark on some issues regarding photometric 
redshift determinations. In order to implement our procedure 
we require more redshift bins than is usual for tomographic 
techniques (which frequently saturate at 2 or 3 bins). This 
requires more accurate photo-z determinations, but still well 
within the range of possible accuracies. For example, SNAP 
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is designed to achieve accurate photo-z measurement with a 
random error a z = 0.03 (Massey et al. 2003). This would 
allow us to perform tomography with up to 10 redshift bins 
in the range < z < 3. Of course it is still crucial to cal- 
ibrate the photo-z estimates to the percent level, since con- 
straints on the cosmological parameters come from measur- 
ing the lensing efficiency to the percent level. This will re- 
quire a spectroscopic survey to calibrate the photo-z distribu- 
tion, as discussed in Bernstein & Jain (2003). Such a mea- 
surement should also mitigate against another important sys- 
tematic arising from the tails of the photo-z distribution: the 
misidentification as source galaxies of a fraction of galaxies 
which actually lie near the lensing plane leading to a mises- 
timation of the lensing efficiency. Both effects can be con- 
trolled introducing gaps, larger than the photo-z error, in the 



galaxy distribution when computing the cross power spectra 
between neighboring bins. We have found that even gaps 
with Az = 0.2, a typical photo-z error from a five color sur- 
vey, enlarge the errors from cross spectrum tomography with 
5 redshift bins by less than 5%, as expected from our earlier 
discussion. While further study (including realistic photo-z 
errors) is needed, this issue is not significantly different for 
cross-spectra and auto-spectra. 
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